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Precision determination of the CKM angle 7/^3 depends upon constraints on charm mixing amplitudes, mea- 
surements of doubly-Cabibbo suppressed amplitudes and relative phases, and studies of charm Dalitz plots 
tagged by flavor or CP cigenstates. In this note we describe the technique used at CLEO-c to constrain the 
Ks,l'K + t{~ model uncertainty, and its impact on 7/03 measurements at S-factories presented at the Charm 2007 
Workshop. 



1. Introduction 



Measurement of the CKM angle 7/03 is challeng- 
ing. Several methods have been proposed using B T — > 
DK T decays; 1) the Gronau-London-Wyler (GLW) 
method [1] where the D decays to CP eigenstates 2) 
the Atwood-Dunietz-Soni (ADS) method [2] where the 
D decays to flavor eigenstates and 3) the Dalitz plot 
method [3, 5] where the D decays to a three-body fi- 
nal state. This latter method has been used recently 
by CLEO to measure the K*K strong phase via the 
three-body decay D° — > K Jr K~ir° [6]. Uncertainties 
due to charm contribute to each of these methods. 
The CLEO-c physics program includes a variety of 
charm measurements that will improve the determina- 
tion of 7/^3 from the i?-factory experiments, BaBar 
and Belle. The pertinent components of this program 
are improved constraints on charm mixing amplitudes 
- important for GLW, measurement of the relative 
strong phase 5 between D° and D° decay to K + tt~ - 
important for ADS, and studies of charm Dalitz plots 
tagged by hadronic flavor or CP eigenstates. The to- 
tal number of charm mesons accumulated at CLEO-c 
will be much smaller than the samples already accu- 
mulated by the i?-factories. However, quantum cor- 
relations in the DD system from -0(3770) provides a 
unique laboratory in which to study charm. 

The decay with the largest branching fraction rel- 
evant to the determination of 7/03 D° — > A's7r + 7r~. 
Recently Babar [7] and Belle [8] have reported 7 = 
(92 ±41 ±11 ±12)° and 3 = (53I^±3±9)°, respec- 
tively, where the third error is the systematic error 
due to modeling of the Dalitz plot. 

Both D° and D° populate the Dalitz plots 



As7r + 7r , (as well 



-tt°, K+K' 



KgK ± ir T ) and so can be used in the determination of 
7/03 which exploit the interference between b — ► cus 
(B~ -> D°K~) and b -> ucs (B~ -> D°K~) where 
the former process is real and the latter is proportional 
to ~ e~ %1 [9]. Studying CP tagged Dalitz plots allows 
a model independent determination of the relative D° 
and D° phase across the Dalitz plot. We describe this 
technique in the following sections. 



2. Determining 7/^3 From B Decays 

Our analysis follows the work outlined in [3], [5], 
and [4]. We consider the decay process B ± — > _DA ± , 
followed by the three-body decay D — ► Kstt + tt" . As- 
suming no CP violation, we define the decay ampli- 
tudes for the D° and D° to be 



A{D° 



;x,y) 



fo(x,y) 
fD(y,x). 



(1) 



Sensitivity to the angle 03 comes from the interference 
of the neutral D mesons from B^ — > DK ± . Since the 
D meson is in a linear combination of flavor states, 
the amplitude for a D° — > Kstt + tt~ event originating 
from a B decay is then 

A{B-^(K S 7T+7T-) D K-) cx f D +r B e ie -f w (2) 

up to an overall normalization. The angle 9 is de- 
fined as 9± = &b ± 03- Here 8b is the strong phase 
difference between color-suppressed and favored am- 
plitudes, whilst tb is the ratio between the color- 
suppressed to favored amplitudes. Theoretical esti- 
mates place tb between 0.1-0.2 [11]. This has been 
confirmed by BaBar (r B = 0.12 ± 0.08 ± 0.03(syst) 
± 0.04(model), [7]) and Belle (r B = 0.16 ± 0.05 ± 
O.Ol(syst) ± 0.05(model), [8]). 

The D° -> K s ir+ir- Dalitz plot is divided into 2Af 
bins, symmetric under exchange of x and y. The bins 
are indexed from — i to i, excluding zero, as in shown 
in Fig 1. The coordinate exchange x <-> y thus corre- 
sponds to the exchange of bins i <-> —i. For simplicity 
we ignore the effects of efficiency and background in 
the Dalitz plot. The number of events in the i-th bin 
of the Kstt + tt~ Dalitz plot from a D decay is then 



and expressed as 



Ki = A D / \f D (x, y)\'dxdy = A D F % . (3) 

The interference between the D° and D° amplitudes 
is parametrized by the two quantities 



Re [f D (x, y)fh(y, x)] dx dy (4) 
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Figure 1: Binning of the D — > Ksir vr Dalitz plot, 
and 

Si = —T===l Im{f D (x,y)fn(y,x)]dx dy, (5) 

V Fit 1 —i JT>i 

where the integral is performed over a single bin. The 
number of events in the i-th bin of the Ksit + it~ Dalitz 
plot from a B decay is then 

N, = K t + r 2 B K- t 

-2r B ^J KiK-i (cj cos 6- - Si sin 6 J) , (6) 

again up to an overall normalization. It is impor- 
tant to note that Ci and Si depend only on the D 
decay. These are the quantities that we measure us- 
ing CLEO-c data. Although in principle they could 
be left as free parameters in a D — > A"s7r + 7r~ Dalitz 
plot analysis from B decays, their values can be more 
precisely determined from correlated Dcp decays pro- 
duced at CLEO-c. 

Thus, we can constrain 0±, and in turn 7/(^3, if we 
know Ki , Ci, and Si . The Ki can be easily determined 
using flavor-tagged D° — > Kstt + it~ Dalitz plot. In the 
next section we show how the Ci can be obtained using 
binned, CP-tagged D° — * Kstt + tt~ Dalitz plots. 

3. Measuring c; From CP-Tagged D 
Decays 

For D mesons that decay into a CP eigenstate, we 
write the initial state of the Dasa linear combination 
of flavor eigenstates 

fcp±(x,y) = -^= \ fo(x,y) ± fo{y,x)\ . (7) 



In terms of this amplitude the number of events in the 
i-th bin of a CP-tagged Dalitz plot is 

Mf = h CP ± (Ki ± 2c, y/KK~ + K^, (8) 

where hcp± is a normalization factor. 

The expression given above for M.f~ can be used to 
measure Ci directly, even if only one type of CP tag 
is reconstructed. Care must be taken to use the cor- 
responding value of hcp± as defined above. However, 
if samples of both CP parities are available we can 
combine the expressions for and M~ to get the 
following equation 

_l (Mr-M+) (Kj + K-j) 
1 2(M? + Mr) JK-R- ' W 

We thus have an expression for measuring a simply by 
counting events within the bins of flavor-tagged and 
CP-tagged Dalitz plots. 

At CLEO-c we produce D°W pairs from the de- 
cay of a tp(3770) in a definite eigenstate of C = — 1. 
Ignoring both the effects of CP violation, the double 
tag rate for final states |1) and |2) is given by 

r(l, 2) = \A(1, 2)| 2 + (mixing terms), (10) 

where 

A(l,2) = (1|P»°)(2|P^) - (l\W)(2\D°) (11) 

For the time being we ignore the effects of correlations 
and mixing in the Kit tagged Dalitz plot. This is not 
expected to make a significant difference for the Kit 
mode, as terms proportional to rK-n — 0.06 and r\^ 
are negligible. 

3.1. Optimized Binning 

Although the quantity Si can only be measured us- 
ing a A's7r + 7r~ vs. Ksit + tt~ double Dalitz analysis [4], 
it can still be approximated from a single Dalitz plot if 
the binning is fine enough. If the bins are small enough 
that the phase difference and the amplitude remains 
constant across each bin, the strong phase parame- 
ters become a = cos(<5d), Si = sin((5£i), so that the 
equality Si = \Jl — (3 is true. It has been shown [3] 
that this equality holds for 200 or more bins, which is 
clearly not feasible for the number of Dcp tags pro- 
duced at CLEOc. In order to circumvent this problem, 
Bondar has proposed an alternate, model-dependent 
method for binning the As7r + 7r~ Dalitz plot [4]. The 
optimal choice depends on the D° — » Ksit + tt~ model. 
In this analysis we use the isobar model amplitude ob- 
tained from the most recent Belle ^3 Dalitz analysis 
[8]. 

From the consideration above it is clear that a good 
approximation to the optimal binning is the one ob- 
tained from the uniform division of the strong phase 
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difference Srj- We thus take the definition of i-th bin 
to be 

2ir(i-l/2)/N < 8 D (x,y) < 2n(i+ l/2)/N. (12) 

An example of such a binning with TV = 8 is shown in 
Fig. 2. 




\L l l l l I l l l l I l l l l 

12 3 

M 2 (K^ + ) 



Figure 2: Divisions of the D° — » Kstt + tt Dalitz plot with 
uniform binning of ASd strong phase difference with N = 
8. 



4. Event Selection 

4.1. Double-Tagged D K s tt + tt~ 
Events 

This analysis uses a combination of two-body CP 
and flavor tags. Since the neutral D mesons are pro- 
duced at ^'(3770) threshold they are correlated in a 
C = — 1 state. If mixing is ignored we can determine 
whether the parent particle was a D° or D°, up to 
DCS contributions. Similarly, if CP violation is ig- 
nored, then the D mesons must be in eigenstates of 
opposite CP [10]. 

To determine the flavor of the D meson, we tag 
D° — > A's7r + 7r~ events with the two-body P° — > 
K + tt~ mode. 1 We use the two CP-even tags K + K~ 
and 7r + 7r~, and the two CP-odd tags Kstt and KsV- 

We introduce two quantities that are reconstructed 
on both sides of a double-tagged decay. T he beam- 
constrained mass is defined as M\, c = y/E^ — p 2 D , 



1 The inclusion of charge-conjugate modes is implied 
throughout our analysis. 



where P& is the beam energy and p D is the square 
of the reconstructed 3-momentum of the D meson. 
We require that the beam-constrained mass of the 
reconstructed candidate is within 3cr of the nominal 
D mass, which corresponds to a selection criteria of 
1.8603 < M bc < 1.8687 GcV. The other quantity is 
the energy difference between the beam and the re- 
constructed D, defined as AE = Eb eam — Ed- We 
apply a selection criteria of |AP| < 30 MeV to all 
D° — * Kstt + tt~ candidates. 

Additional selection criteria are placed on the du- 
aghter particles to ensure basic track quality. For ex- 
ample, we select pion track momenta between 0.05 < 
p < 2.0 GeV. Both signal and tagging modes contain- 
ing a Ks are selected to be within 3<r of the Ks mass, 
which corresponds to ±7.5 MeV from the central Ks 
mass value of 497.6 MeV. 

We only reconstruct Ks particles that decay 
through the n + ir~ channel; we do not attempt to re- 
construct Ks — > 7r 7r°. Fake Ks candidates can be 
misrcconstructcd from combinatoric 7r + 7r~ pairs. To 
suppress these events we apply a selection criteria on 
the flight significance f s > to our Ks candidates. 
Additionally, we require that the it mass falls within 
3cr of its nominal value. 



4.2. Double-Tagged D — > K l tt+7t- 
Events 

For D° — > Al7t + 7t~ decays we require the same 
selection criteria on charged pions and ir° candidates 
as those described for D° — > Ksir + n~ decays. How- 
ever, because of the large flight distance of the K L , 
the A£7r + 7r~ signal is reconstructed using a missing 
mass technique. We require the signal side to have 
exactly two charged tracks. We also apply ir° , 77, and 
Ks vetoes. Using the measured momentum of the 
tagged D, we compute the missing momentum and 
energy on the signal side. We require that the missing 
mass squared satisfies the condition 0.21 < m? < 0.29 
GeV 2 . The background for D° — > Klit+tt^ mode is 
approximately 5%. 



4.3. Double-Tagged K L ir () vs. K s tt + tt~ 

We can increase our statistics by reconstructing 
D a — > A's7T + 7r _ events tagged with the CP-even 
mode K^iT . Wc require zero tracks and exactly one 
7T° candidate on the tag side. We veto events contain- 
ing 77 candidates, and impose similar criteria on the 
Kl missing mass as described above. 

The final yields for all tag modes are summarized 
in Table I 
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Table I Yields for CP-tagged K s tv + it~ and A'l7t + 7t in 
398 pb -1 data, by tag mode. 



Tag Mode 






K+K~ 


61 


194 


7T + 7T~ 


33 


90 




108 


263 


Ksn 


29 


21 


K l -k° 


190 





5. Combining AVr + 7r and K L ^ 



7T 



The tagged Klit+tt^ Dalitz plots are included to in- 
crease the statistical accurancy of this analysis. How- 
ever, if we naively combine the Dalitz plots with K$ 
and Kl we will find our measurement of Cj to be bi- 
ased. We must first account for the phenomenological 
differences between the Kstt + t:~ and A'L7r+7r~ mod- 
els. 

Since the K$ and Kl mesons are of opposite CP, 
the doubly-Cabibbo suppressed amplitudes in each 
Dalitz plot will contribute with opposite signs. We 
can see this by inspecting the D° decay amplitude for 
each each Dalitz plot 



A(Ksttit) 
A(K l tttt) 



1 
1 

71 



A{K°tttt) + A{K°mr) 

A(K%Tr) - ACRQinr)] (13) 
(14) 



The effect of this relative minus sign is to intro- 
duce a 180° phase for all DCS K* resonances in the 
K model. We can use [/-spin symmetry to 
relate the amplitudes for resonances of definite CP 
eigenvalue. We find that these states aquirc a factor 
of rxe lSl< — — tan 0c, where 9c is the Cabibbo angle. 

In our study we mulitply all DCS amplitudes in the 
Kltt + tt~ model by -1. From this "base" model we fix 
rx = 0.06 for each CP eigenstate, then vary the phase 
8k between and 2tt. For each bin we then find the 
largest resulting deviation in ct, and report this value 
as the systematic uncertainty in the KLir + ir~ model. 

To better understand the difference between the 
A'57r + 7r~ and Kltt^~tt~ models, we compare the nu- 
merically calculated values of Cj in each Dalitz plot. 
Wc find that the value for Ci is systematically larger in 
each bin for Kltt + tt^. In Fig. 3 we can see that the 
difference is significantly larger than the systematic 
uncertainty in our ATi,7r + 7r~ model. 




,1 Ll- t I 1 I J L I J I J I i I I J I I 1 I I L I I L I L I I I I i I I I I J I 

1 2 3 4 5 6 7 8 

Bin Number 



Figure 3: Values for d numerically determined from our 
model. Az,7r + 7r - values are in red, Kstt + -k~ in blue. In 
each bin the values for a are systematically larger for 
Kltv + tv~. The error bars represent the uncertainties in 
the Kltt + tt~ model parametrization. 



6. Results 

We report the difference in Ci between Kstt + tt~ 
and Kt,k + tt~ as measured in 398 pb" 1 of data. In 
Fig. 4 we compare the Ci differences calculated from 
our model and measured from data. The error bars 



0.8- 




g 7 a 

Bin Number 



Figure 4: The difference in Ksirir and Kl^tt values of a, 
numerically determined from our model (black) and mea- 
sured in 398 pb" 1 of data (green). The green error bars 
represent the combined statistical and model uncertainty. 



in this figure represent both statistical and model un- 
certainty combined. With a reasonable understanding 
of the Ci between the Kstt + tt~ and Kltt + tt^ Dalitz 
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plots, we can estimate the final precision with which 
we expect to measure c, once 750 pb" 1 of data is avail- 
able. The values of Cj from our study are once again 
plotted in Fig. 5, but here the error bars represent the 
statistical uncertainty obtained from 398 pb _1 of data 
scaled up to 750 pb _1 . 
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Figure 5: The central values of a are computed from our 
model with expected sensitivity from 750 pb -1 of data. 
The error bars are determined by scaling the statistical 
uncertainty obtained from 398 pb _1 of data, then combin- 
ing the Kltt + tv~ model uncertainty. 

We expect good sensitivity to the measurement of Ci 
with the entire CLEO-c data. This measurement can 
reduce the model uncertainty on 7/(^3 to a precision 
of about 4° [3]. 
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